LGP-30 USERS' ORGANIZATION - POOL 

IrArod^ic jygiy Information for Subrouti nes No. F2-38a an d F2-38b. 

These two related but operatively independent subroutines are useful 
for two types of jointly distributed variables, depending on their mean values. 
Subroutine No. F2-38a is appropriate if either mean value is non-zero, and the 
generated least squares regression line has the form, 

I c = a + bX 

The point defined by the mean values of the variables, (X,Y)lies on the line, 
either (or both) of the mean values are non-zero, and the line does not pass 
through the origin (or what is the same thing, the constant "a" is non-zero). 
If, on the other hand, the bivariate distribution's variables happen to 
have mean values which are both zero or, if a transformation is made on X and Y 
such that the transformed variables, x and y, meet this condition, then Sub- 
routine No. F2-38b is fitting. It may be used to fit the regression line 

y c - b x 

to the joint distribution. 

The following symbol-conventions are used: 

1. Capital alphabetical symbols, X, and Y, represent "observed" 

data or "scores" measured from an untampered origin; i.e., no trans- 
formation of the type listed below has been made on such X, Y data. 

2. Small alphabetical symbols, x, and y, are used to represent deviations 
or measurements from sample mean-values X", and Y, and may be obtained 
from initial X, and Y measurements or scores by the relations: 

x = X - X 

y = Y - Y 

3. The Y (or y) symbols with subscripts "c" represent "estimated" or 
computed values supplied by the least squares equations. 

Both No. F2-38a and F2-38b use the original floating point interpretive 
routine, No. 2lj..O and its input-output companion, No. 11.6-12.6, 

These regression programs come equipped with tapes having three short 
examples and a control program for directing the computer through its operations 
while executing the computations. Information is punched in the following order 
with one tape each for No. F2-38a, and No. F2-38b: 

1. Bivariate Linear Regression (Either F2-38a or F2-38b) 

2. Short Control Program 

3. Three examples 

The Flexowriter printing done while the computer executes these sets of examples 
appears on pages k and 12 of this program description. In practice, regression- 
correlation calculations would not be made on such small sets of data, but these 
serve well to illustrate program operation. It will also be noted that the 
usual equal-spacing requirement for simplified data-transformation has not been 
met | this condition was down-graded when an interest developed in applying 
Mo. F2-38b to No, F2-38a's previously constructed examples. 

c— 
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Program No. F2-38a 
Dec. 1, 1959. 

Program Title: "Bivariate Linear Regression" 

(Straight Line Fit to Observed Data) 

Author: R. L. Stearman 

Installation: Booz-Allen Applied Research, Inc., Bethesda, Maryland. 

Purpose: To compute the least squares regression line, of the form 

Y Q = a + bX, to a set of points in floating point and to compute 
the correlation coefficient, r, for this set of points. (See 
list of additional statistics computed, as given below under 
"Output"). 

Input: (Note: This is the same for both F2-38a and F2-38b) 

1. Location of the floating point interpretive 
subroutine, No. 2i|.0. 

2. Location of N, the size of sample (N must be given 
in floating point). 

3. Location of the first value of X (values of X must 
be in floating point and stored sequentially) . 

k* Location of the first value of Y (values of Y must 
also be in floating point and stored sequentially 
according to the order of the corresponding values of X). 

Subroutine must be entered in fixed point and exits in fixed point. 

Output : 

The following statistics are stored within the program: 

Statistic Storage Location 

JC X L + 0103 



o 



EX 2 L Q + 0101; 

£x » s .s. x =J> - gx) 2 L q + 0106 

L Y L + 0107 

o 

S Y L + 0108 



T ~ ~y =C? - (YXi -u 



Zu y = s.s. v =/Jy _ Q_j) Ln + 0110 

N 



^XY L + 0111 

E xy = s.p.^ =T;XY - (Ex) (£y) L + 0113 

*> * S-P'xy/ S.S. x , L + Olli; 

a = £ T " bgx) L + 0115 

— r— (£) 
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Output; Continued 

Statistic * ( See definitions below) 



Program No. F2-38a 



Storage Location 



r = s.p 



•*/y 



s • s* s# s« 

x y 



z i 

S.S. reg = S «P'xy/ S,S *X 



s.s. = S .S. ~ S » S « 

y* x y reg 



d.f. a n 

2 
s 

y.x 



s.s. y#x / d.f. 



F = s.s. / s % r „ 
reg' y.x 

Lin kage & Calling Sequence: (Note: this is the 

F2-38a and F2-38b) 

To be entered while in fixed point. 



L o 


+ 


0102 


L 




+ 


0116 


L o 


+ 


0117 


L 




+ 


0118 


L 




+ 


0119 


L 




+ 


0120 



same for both Subroutines 



Location Instruction 


Address 


-t • - - • R ..-- 


- L 




0^+1 U - - - - 


L 
... 


iV- +2 . _ Z 


Location of floating 
point 


06+3 - z 


Location of n 


c* + i Z 


Location of first x 


<X + 5 z 


Location of first y 


Cy. + 6 etc. 




Storage: 




2 tracks and 2li sectors 


01 I 


Program Stops: 


None 




-x- Further Definition of Expressions. 





r : Correlation Coefficient. 

s.s. : Sum of Squares due to regression 

s.s. x : Sum of Squares due to residual 

d.f. : Degrees of Freedom 



y.x 



Mean Square due to^residual 



o 
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"Bivariate Linear Regression" Program No. F2-38a 
(Straight Line Fit to Observed Data) 

CONTROL PROGRAM FOR OPERATING No. 38a and THREE SHORT EXAMPLES 



CONTROL 
PROGRAM 



4 



/I0003 too • /ooo3 I *oo ' 



EXAMPLE 

NO. I 



xrl200' 
xul200' 
xiOOOO 1 
xeOOOO' 
xr3000' 
xu3000« 
xzl200' 
xz^OOO' 
xzUooV 
xzU032' 
xzOOOO ' 
xpl600' 
xzOOOO' 
xrl200' 
xul200' 

xb3115 f 
xpOOOO ' 
xb311 J 4-' 
xpOOOO ' 
xb3l02' 
xpOOOO ' 
xb3120« 
xpOOOO ' 
, xeOOOO' 
\xzOOOO* 

^.0003^00 • 
+oo4ooo f 
o'i^^'S 1 

+00^)32' 



5' 



lo. floating point 

input data 

exit f.p. 

L-zero Bivariate Linear Regression 

floating point lo. 

lo. n 

lo. X s 

lo. Y s 

stop 

c.r. 

delay 



a 

print a 
b 

print b 
r 

print r 
F 

print F 
exit f.p. 
stop 



X s 



-0000000 • 



Y s -0000000' 



\ .i6k6oiQ oi .3761061 00 .7996127 00 . 533.9016 01 



. 0003^00 • 

+00^000' 

5' 
o'l^^'V 

EXAMPLE < + 0lli032 ! 



NO. II 



0'i9 , Uo , 59 , 8l' 
V .6000023- 01- 



n 

X s -0000000' 

Y s -0000000' ' 

.2020000 01 .9997060 



00 



.5101681 Ok 



EXAMPLE 
NO. Ill 



/ . 0003*400' 
+001*000' 6' 
52 » 110 ' 205 ' 377 ' Uoo ' 500 ' 
+00U032' 
252 * 280 ' 360 ' 1*00 « if-50 « 520 » 



-0000000' 
-0000000' 



+0000000' 



V^.2236380 03 .5597153 00^ .98069^0 00 



00 



.1006043 03 
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Program No. F2-38a 
Dec. 1, 1959. 

Author; R„ L. Stearman 

Installation: Booz-Allen Applied Research Inc., Bethesda Ik, Maryland. 



Entry- 



Floating Point 
"BIVARIATE LINEAR REGRESSION" 
(Straight Line Fit to Observed Data) 



J. 



Initialize 
transfer to 
floating point 



Counter i ^S 



<? j finished?!' 



f 



-.-vZ&~r, 



Initialize 

program 
return 



i 



Transfer calling | 
sequence 
locations to 
loop initial! zine, 
orders 



* 



NO 



Add X. to 

yx. 



* 



Add X, to 



Initialize loop 
orders for 



~4 



, I Add Y. to 



.1 



± 



Bring N to 
the counter! 



|£ T i 



4<- 



Add Y, to 



± 



JL 



Clear summation 
•storage to zero 



Add X ± Y. to 



JsL 



Transfer to 
floating point 



I 




V 



Modify 
commands 
in loop 



Compute and 



store ss 



X 



i.l 



Compute and 
store ss 

y 



v- 



Compute and 
store sp 



xy 



T* 



! Compute and! 
store b ,• 

i 
V 



Compute and 
store a 



* 



Compute and 
store d. f. 



i 



Compute and 

store ss 

reg 



± 



Compute and 
store ss 

y.x 



V 



Compute and 

store s 

y.x 



Return to 
program 



< 



1 ^™ 

[Compute and' 
store F 









J 




Exit floating 
point 


<— 


V ... 




Compute and 
store r 







© 
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